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A  method  for  modelling  the  thermal  behaviour  of  a  lithium-ion  battery  (LIB)  during  charge  is  presented. 
The  effect  of  charge  conditions  on  the  thermal  behaviour  is  examined  by  means  of  the  finite  element 
method.  A  comparison  of  the  experimental  charge  curves  with  the  modelling  results  validates  the  two- 
dimensional  modelling  of  the  potential  and  current  density  distribution  on  the  electrodes  of  an  LIB  as  a 
function  of  charge  time  during  constant-current  charge  followed  by  constant-voltage  charge.  The  heat 
generation  rates  as  a  function  of  the  charge  time  and  the  position  on  the  electrodes  are  calculated  to 
predict  the  temperature  distributions  of  the  LIB  based  on  the  modelling  results  for  potential  and  current 
density  distributions.  The  temperature  distributions  obtained  from  the  modelling  are  in  good  agreement 
with  the  experimental  measurements. 
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1.  Introduction 

The  lithium-ion  battery  (LIB)  is  the  preferred  power  source  for 
hybrid  electric  vehicles  (HEVs)  and  electric  vehicles  (EVs)  due  to 
its  high  specific  energy,  high  voltage,  and  low  self-discharge  rate. 
To  provide  affordable  HEVs  and  EVs  with  satisfactory  performance, 
there  are  still  major  technical  challenges  to  improve  the  charac¬ 
teristics  of  LIBs  such  as  low  cost,  high  specific  energy,  high  power 
density,  long  service  life,  and  proven  safety.  Because  safety  is  one 
of  the  major  challenges  for  the  LIB  in  HEV  and  EV  applications,  the 
thermal  control  of  LIBs  is  especially  important.  The  main  concern 
in  the  thermal  control  of  LIB  is  the  significant  temperature  increase 
that  can  occur  during  high  power  extraction  and  rapid  charging  in 
HEV  and  EV  applications,  which  may  cause  battery  degradation  and 
thermal  runaway.  Thermal  modelling  can  play  a  key  role  in  control¬ 
ling  the  operating  temperature  and  temperature  uniformity  of  LIBs 
within  a  suitable  range  [1-3]. 

There  have  been  many  previous  studies  of  the  thermal  mod¬ 
elling  of  LIBs  [4-21  ].  Bernardi  et  al.  [4]  developed  a  general  energy 
balance  for  battery  systems  including  the  contributions  from  elec¬ 
trochemical  reactions,  mixing  enthalpies  and  phase  changes.  Rao 
and  Newman  [5]  presented  a  method  of  calculating  the  heat  gen¬ 
eration  rate  based  on  the  general  energy  balance  and  the  enthalpy 
potential  method  for  a  LIB  system.  Chen  and  Evans  [6-8]  developed 
two  and  three  dimensional  models  to  study  the  thermal  behaviour 
of  lithium  polymer  batteries  (LPBs)  and  LIBs.  They  assumed  that 
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the  heat  generation  rate  is  uniform  throughout  the  cell.  Pals  and 
Newman  presented  a  one-dimensional  model  for  predicting  the 
thermal  behaviour  of  LPBs  for  a  single  cell  [9]  and  a  cell  stack  [10]. 
Verbrugge  [11]  modelled  the  three-dimensional  current  and  tem¬ 
perature  distributions  in  LPB  modules.  Botte  et  al.  [12]  studied  the 
influence  of  design  variables  on  the  thermal  behaviour  of  LIBs  based 
on  a  one-dimensional  model.  Al-Hallaj  et  al.  [13]  presented  simpli¬ 
fied  one-dimensional  thermal  modelling  with  lumped  parameters 
to  simulate  the  temperature  profiles  inside  LIB  cells.  Song  and  Evans 
[14]  developed  an  electrochemical-thermal  model  of  LPBs  by  cou¬ 
pling  a  two-dimensional  thermal  model  with  a  one-dimensional 
electrochemical  model.  Gu  and  Wang  [15]  and  Srinivasan  and 
Wang  [16]  developed  a  two-dimensional  thermal  and  electrochem¬ 
ical  coupled  model  to  analyze  the  electrochemical  and  thermal 
behaviour  of  lithium-ion  battery  cells.  Gomadam  et  al.  [  1 7]  reduced 
a  two-dimensional  heat  conduction  model  to  a  one-dimensional 
model  for  LIBs  with  a  spiral  geometry  by  a  coordinate  transfor¬ 
mation  approach.  Chen  et  al.  [18]  developed  a  three-dimensional 
model  of  LIBs  considering  the  location-dependent  convection  and 
radiation  at  the  boundaries  to  reflect  the  different  heat  dissipation 
performances  on  all  of  the  surfaces.  Kim  et  al.  [19,20]  presented 
two-dimensional  thermal  modelling  based  on  the  potential  and 
current  density  distributions  obtained  by  the  same  procedure  as 
that  described  by  Kwon  et  al.  [21  ]. 

Most  of  the  thermal  modelling  studies  of  LIBs  and  LPBs  men¬ 
tioned  above  [4-6,9]  focused  on  the  thermal  behaviour  during 
discharge  and  only  a  couple  [7,8]  briefly  dealt  with  the  thermal 
behaviour  during  charge  to  accommodate  the  dynamic  cycling  of 
batteries  for  EV  applications.  Although  the  heat  generation  rate 
during  charge  is  less  than  that  during  discharge  at  the  same  rated 
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current  of  the  battery,  the  increase  in  local  temperature  becomes 
significant  during  rapid  charge  especially  for  the  large-size  LIBs 
used  in  EV  and  HEV  applications  [3].  Because  there  can  be  many 
circumstances  under  which  rapid  charge  becomes  a  necessity  for 
EV  and  plug-in  HEV  applications  [22-24],  it  is  important  to  develop 
a  reliable  modelling  methodology  to  predict  accurately  the  thermal 
behaviour  of  LIBs  under  various  charge  conditions  for  the  purpose 
of  improving  their  performance  and  life,  as  well  as  ensuring  ther¬ 
mal  safety.  In  this  work,  two-dimensional  modelling  is  performed 
to  calculate  the  potential  and  current  density  distributions  during 
charge  on  the  electrodes  of  an  LIB  using  a  procedure  similar  to  that 
reported  by  Kwon  et  al.  [21  ].  The  modelling  is  validated  by  compar¬ 
ing  the  experimental  and  modelled  charge  curves  at  charge  rates  of 
1,  3,  and  5  C.  Then,  the  thermal  behaviour  is  modelled  as  a  combi¬ 
nation  of  the  heat  generation  due  to  the  electrochemical  reactions 
and  electrical  resistance  within  the  cell  as  a  function  of  the  electric 
current  and  state-of-charge  (SoC)  using  a  similar  procedure  to  that 
described  by  Kim  et  al.  [19,20].  The  electrical  heating  rate  is  cal¬ 
culated  from  the  electrical  current  flow,  which  can  be  evaluated 
based  on  the  potential  distributions  on  the  electrodes  obtained 
from  the  electrochemical  modelling.  The  thermal  model  is  vali¬ 
dated  by  comparing  the  modelling  results  with  the  experimental 
measurements. 


Table  1 

Fitting  parameters  used  to  calculate  potential  and  current  density  distributions  on 
electrodes. 


Parameter 

Value 

a0  (V) 

4.135 

at  (V) 

-0.91 

a2  (V) 

1.405 

a3  (V) 

-1.48 

a4  (Am-2) 

116.859 

a5  (Am-2) 

-892.8001 

a6  (Am-2) 

5250.46 

07  (Am-2) 

-13,623.09 

a8  (Am-2) 

15,853.17 

ag  (Am-2) 

-6757.8539 

Table  2 

Parameters  used  for  thermal  modelling. 


Component 

Density 

Heat  capacity 

Thermal 

(gem-3) 

(Jg-’K-1) 

conductivity 
(Wcm-1  K-1) 

Current-collector  of 

2.7 

0.9 

2.38 

positive  electrode 
Electrode  material  of 

1.5 

0.7 

0.05 

positive  electrode 
Current-collector  of 

8.96 

0.385 

3.98 

negative  electrode 
Electrode  material  of 

2.5 

0.7 

0.05 

negative  electrode 
Separator 

1.2 

0.7 

0.01 

Pouch 

1.15 

1.9 

0.16  x  10-2 

2.  Mathematical  model 

A  14.6  Ah  LIB  comprising  a  LiMn204  cathode,  a  graphite  anode, 
and  a  plasticized  electrolyte  from  LG  Chem.  was  modelled  in  this 
work.  A  cell  consisting  of  two  parallel-plate  electrodes  of  the  bat¬ 
tery  shown  in  Fig.  1  was  chosen,  because  the  battery  consisted  of 
the  same  repeating  units  of  positive  and  negative  electrode  plates, 
polymer  electrolytes  and  separators.  The  tabs  are  the  current- 
collectors  extending  outside  from  the  rectangular  electrodes  and 
they  do  not  contain  the  electrode  (active)  material.  A  schematic  dia¬ 
gram  of  the  current  flow  in  the  cell  during  charge  is  also  included 
in  Fig.  1.  The  distance  between  the  electrodes  was  assumed  to  be 
so  small  that  the  current  flow  between  the  electrodes  would  be 
perpendicular  to  them.  The  modelling  procedure  used  to  calculate 
the  potential  and  current  density  distribution  on  the  electrodes 
was  similar  to  that  used  by  Kwon  et  al.  [21].  From  the  continu¬ 
ity  of  the  current  on  the  electrodes,  the  following  equations  could 
be  derived: 
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Fig.  2.  Dimensions  of  electrodes  and  positions  of  tabs  of  a  1 4.6  Ah  LIB  from  LG  Chem. 
Dimensions  in  mm.  Arrow  passing  through  tab  from  positive  electrode  represents 
a  vertical  line  along  which  temperature  distributions  are  plotted  in  Fig.  9. 
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Fig.  3.  Schematic  illustration  of  discharge-charge  protocol  of  battery. 

V  •  ip  -\-J  =  0  in  £2p  (1) 

V  •  in  —  J  —  0  in  £2n  (2) 

where:  ip  and  in  are  the  linear  current  density  vectors  (current 
per  unit  length  (Am-1))  in  the  positive  and  negative  electrodes, 
respectively;  J  is  the  current  density  (current  per  unit  area  (A  m-2 )) 
transferred  through  the  separator  from  the  negative  electrode  to 
the  positive  electrode;  £2P  and  £2n  denote  the  domains  of  the  posi¬ 
tive  and  negative  electrodes,  respectively.  According  to  Ohm’s  law, 
ip  and  in  can  be  written  as: 

ip  =  -  — VVp  in  C2p  (3) 

rP 

in  —  —  ~  in  £2n  (4) 

Pn 

where:  rp  and  rn  are  the  resistances  {Q)  of  the  positive  and  nega¬ 
tive  electrodes,  respectively;  Vp  and  Vn  are  the  potentials  (V)  of  the 
positive  and  negative  electrodes,  respectively.  The  resistances,  rp 
and  rn,  were  calculated  as  described  in  Refs.  [19-21  ].  The  following 
Poisson  equations  for  Vp  and  Vn  were  obtained  by  substituting  Eqs. 
(3)  and  (4)  into  Eqs.  (1)  and  (2),  respectively: 


V2Vp  —  J-Ppj  in 

(5) 

V2Vri  =  —PnJ  in  £2n 

The  relevant  boundary  conditions  for  Vp  were: 

(6) 

dVD 

lif-0  r« 

(7) 

on  r„ 

Vp  dn  L 

(8) 

where:  9/9  n  denotes  the  gradient  in  the  direction  of  the  outward 
normal  to  the  boundary.  The  first  boundary  condition  (7)  implies 
that  there  is  no  current  flow  through  the  boundary  ( rp  i )  of  the  elec¬ 
trode  other  than  the  tab.  The  second  boundary  condition  (8)  means 
that  the  linear  current  density  through  the  tab  (rp2)  of  length  L 
(cm)  is  constant  with  a  value  of  /0/L,  where  /0  is  the  total  current  (A) 
through  the  tab  in  the  mode  of  constant-current  charge,  for  which 
the  boundary  condition  at  rp2  becomes  simply  Vp  =  Vc  where  Vc  is 
the  fixed  value  of  charge  voltage.  The  boundary  conditions  for  Vn 


are: 

TiH  “  r". 

0) 

Vn  =  0  on  r„2 

(10) 

The  first  boundary  condition  (9)  implies  the  same  as  in  the  case  of 
Vp.  The  second  boundary  condition  (10)  means  that  the  potential  at 


the  tab  of  the  negative  electrode  has  a  fixed  value  of  zero  as  the  ref¬ 
erence  potential.  The  difference  between  the  governing  equations 
for  charge  used  in  this  work  and  those  for  discharge  is  that  the  signs 
in  front  of  J  in  Eqs.  (1),  (2),  (5)  and  (6)  are  opposite  to  those  in  the 
corresponding  equations  of  the  previous  studies  [19-21  ].  The  only 
difference  between  the  boundary  conditions  for  charge  used  in  this 
work  and  those  for  discharge  is  that  the  sign  in  front  of  Jo  in  Eq.  (8) 
is  opposite  to  that  in  the  corresponding  boundary  condition  of  the 
previous  investigations  [19-21]. 

The  current  density,  J,  of  Eqs.  (5)  and  (6)  is  a  function  of  the 
potential  difference  between  the  positive  and  the  negative  elec¬ 
trodes  {Vp  -  Vn).  The  functional  form  depends  on  the  polarization 
characteristics  of  the  electrodes.  In  this  work,  the  following  polar¬ 
ization  expression  used  by  Tiedemann  and  Newman  [25]  and 
Newman  and  Tiedemann  [26]  was  adopted 


J  =  Y{VP  -Vn-U)  (11) 

where  Y and  U are  the  fitting  parameters.  As  suggested  by  Gu  [27],  U 
and  y  were  expressed  as  functions  of  the  depth-of-discharge  (DoD) 
as  follows: 

U  =  a0  +  ai(DoD)  +  a2(DoD)2  +  a3(DoD)3  (12) 

Y=a4+a5(DoD)+a6(DoD)2+a7(DoD)3+a8(DoD)4  +  a9(DoD)5  (13) 

where  a0-a9  are  constants  to  be  determined  experimentally.  The 
fitting  parameters  used  to  calculate  the  potential  and  current  den¬ 
sity  distribution  on  the  electrodes  are  listed  in  Table  1. 


Fig.  4.  Variations  of  (a)  charge  voltage  and  (b)  charge  current  during  CC-CV  charg¬ 
ing. 
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0.0014 

0.0013 

0.0012 
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0.0009 

0.0008 

0.0007 
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0.0004 
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Fig.  5.  Distributions  of  (a)  potential  on  positive  electrode,  (b)  potential  on  negative  electrode  and  (c)  current  density  on  electrodes  of  LIB  at  charge  time  of  45  min  during  CC 
charging  at  1  C. 


By  solving  the  equations  listed  previously,  the  distribution  of  the 
current  density,  J,  on  the  electrodes  can  be  obtained  as  a  function 
of  the  position  on  the  electrode  and  time.  Therefore,  the  DoD  varies 
with  the  position  on  the  electrode  and  the  time  elapsed  during 
charge.  The  distribution  of  DoD  on  the  electrode  could  be  calculated 
from  the  distribution  of  J  as  follows: 


DoD  =  DoD0  - 


(14) 


where:  DoD0  is  the  initial  value  of  DoD;  t  is  the  charge  time  (s);  Qj 
is  the  theoretical  capacity  per  unit  area  (Ahm-2)  of  the  electrodes. 

The  thermal  modelling  procedure  was  used  to  calculate  the  tem¬ 
perature  distribution  on  the  electrodes  was  similar  to  that  used  by 
Kim  et  al.  [19,20].  Since  the  thickness  of  the  battery  cell  was  much 
shorter  than  its  other  dimensions,  the  temperature  variation  along 
the  z  direction  in  Fig.  1  could  be  neglected.  Based  on  the  differential 
energy  conservation  for  a  battery,  the  transient  two-dimensional 
equation  of  heat  conduction  was  written  as  follows: 


pCp 


dr 

at 


+  Q  —  Qconv 


(15) 


where:  pis  the  density  (kg  m-3);  Cp  is  the  volume  averaged  specific 
heat  capacity  at  constant  pressure  (Jkg-1  °C-1);  T  is  the  tempera¬ 
ture  (°C);  kx  and  ky  are  the  effective  thermal  conductivities  along 
the  x  and  y  directions  (refer  to  Fig.  1  for  the  x  and  y  directions) 
(Wm-1  °C-1 ),  respectively,  q  is  the  heat  generation  rate  per  unit 
volume  (W  m-3 );  qconv  is  the  heat  dissipation  rate  (W  m-3 )  through 
the  surfaces  of  the  battery  by  convection.  The  effective  thermal 
conductivities  of  the  various  compartments  of  the  cell  could  be 


estimated  based  on  the  equivalent  networks  of  parallel  and  series 
thermal  resistances  of  the  cell  components  [7,18]. 

The  heat  generation  rate,  q ,  was  given  as  follows: 


q  =  aj 


Eoc-E-T 


dE0C 

~df 


+  aprpip  +  anrni„ 


(16) 


where:  a  is  the  specific  area  of  the  battery  (m-1);  J  is  the  current 
density  (Am-2)  calculated  by  Eq.  (8);  Eoc  is  the  open-circuit  poten¬ 
tial  of  the  cell  (V),  E  is  the  cell  voltage  (V),  ap  and  an  are  the  specific 
area  of  the  positive  and  negative  electrodes  (m-1 ),  respectively;  ip 
and  in  are  the  magnitudes  of  the  vectors  ip  and  in  obtained  by  Eqs.  (3 ) 
and  (4)  (Am-1),  respectively.  The  detailed  definition  of  each  term 
on  the  right-hand  side  of  Eq.  (16)  was  described  in  Refs.  [16,19,28]. 
The  heat  dissipation  rate,  qCOnv,  was  derived  as  follows: 

qconv  =  —  Tair)  (17) 

where:  h  is  the  convective  heat- transfer  coefficient  on  the  surfaces 
of  the  battery  (Wm-2oC-1);  d  is  the  thickness  of  the  battery  cell 
in  the  direction  perpendicular  to  the  parallel  electrodes  (m);  Tair  is 
the  ambient  temperature  (°C).  The  rate  was  rendered  by  approx¬ 
imating  a  three-dimensional  object  to  a  two-dimensional  one,  as 
shown  in  Eq.  (15).  The  convective  boundary  condition  applied  to 
the  boundaries  of  the  electrode  is  reported  by  Kim  et  al.  [19].  The 
parameters  used  for  the  thermal  modelling  are  listed  in  Table  2. 


3.  Results  and  discussion 

The  solutions  to  the  governing  Eqs.  (5),  (6)  and  (15)  subject  to 
the  associated  boundary  conditions  were  obtained  using  the  finite 
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Fig.  6.  Temperature  distributions  based  on  experimental  IR  image  and  modelling 
for  LIB  at  charge  times  of  (a)  8  min,  (b)  1 6  min  and  (c)  36  min  during  CC-CV  charging 
with  CC  charge  at  3  C. 


element  method.  The  finite  element  package  used  is  an  in-house 
code  specifically  developed  to  solve  the  problem  defined  in  this 
work.  In  order  to  test  the  validity  of  the  model,  experiments  were 
carried  out  at  a  room  temperature  of  25  °C  using  a  14.6  Ah  LIB  fab¬ 
ricated  by  LG  Chem,  of  which  the  dimensions  of  the  electrodes  and 
the  positions  of  the  tabs  are  shown  in  Fig.  2.  The  experimental  pro¬ 
cedures  are  illustrated  schematically  in  Fig.  3.  To  maintain  the  same 
states-of-charge  before  the  beginning  of  the  charge  experiments, 
the  battery  was  discharged  under  the  constant-current  (CC)  rate  of 
1  C  until  the  voltage  dropped  to  the  cut-off  voltage  of  3.0  V  and  then 
it  was  allowed  to  rest  for  8  h.  The  battery  was  charged  with  CC  rates 
of  1,  3,  and  5  C  until  the  voltage  reached  the  predetermined  value 
of  4.15  V.  Then,  the  charging  voltage  was  kept  constant  at  4.15  V 
until  the  charge  current  tapered  to  0.75 A  and  this  charging  mode 
is  referred  to  herein  as  constant-voltage  (CV)  charging. 

The  variation  of  the  voltage  and  current  of  the  battery  dur¬ 
ing  CC-CV  charging,  is  shown,  in  Fig.  4(a)  and  (b).  The  modelling 
results  and  experimental  measurements  are  in  good  agreement 
with  each  other.  The  model  predicts  a  decrease  in  the  CC  charg¬ 
ing  time  and  an  increase  in  the  CV  charging  time  with  increasing 
rates  of  CC  charging.  The  distributions  of  potential  and  current  den¬ 
sity  on  the  electrodes  during  charge  were  obtained  as  a  function  of 
time  for  various  charge  rates.  As  an  example,  Fig.  5(a)-(c)  shows 
the  distributions  of  the  potential  on  the  positive  electrode,  the 
potential  on  the  negative  electrode,  and  the  current  density  at  a 
charge  time  of  45  min  during  CC  charging  at  1  C,  respectively.  In 
Fig.  5(a),  the  potential  gradient  on  the  positive  electrode  appears 
to  be  most  severe  in  the  region  where  the  tab  is  attached  to  the 
current-collector.  This  is  because  all  of  the  current  flows  through 
the  conducting  current  collector  from  the  tab  into  the  entire  elec¬ 


trode  plate.  Again,  the  potential  gradient  on  the  negative  electrode 
shown  in  Fig.  5(b)  is  the  highest  in  the  region  near  the  tab,  because 
all  of  the  current  needs  to  flow  from  the  entire  electrode  plate  into 
the  tab.  The  non-uniform  distributions  of  current  density  trans¬ 
ferred  from  the  positive  electrode  to  the  negative  electrode  during 
CC  charging  at  1  C  are  presented  in  Fig.  5(c). 

After  obtaining  the  distributions  of  potential  and  current  den¬ 
sity  on  the  electrodes  during  charge,  the  temperature  distributions 
of  the  battery  can  be  calculated  as  a  function  of  time  for  vari¬ 
ous  charge  rates  using  Eq.  (15).  As  a  demonstration,  Fig.  6(a)-(c) 
shows  the  temperature  distributions  based  on  the  experimental  IR 
image  and  the  modelling  for  a  14.6  Ah  LIB  at  charge  times  of  8, 16 
and  36  min  during  CC-CV  charging,  where  the  rate  of  CC  charg¬ 
ing  is  3  C.  Similar  plots  to  Fig.  6(a)-(c)  are  made  in  Fig.  7(a)-(c) 
for  charge  times  of  3,  6  and  24  min,  respectively,  during  CC-CV 
charging,  where  the  rate  of  CC  charging  is  5  C.  In  the  IR  images, 
the  region  within  the  yellow  lines  of  the  rectangle  corresponds  to 
the  computational  domain  for  the  modelling.  Therefore,  the  max¬ 
imum  temperature  within  the  computational  domain  should  be 
differentiated  from  that  of  the  current  collecting  tab  protruding 
beyond  the  electrode  region  of  the  battery  cell.  The  overall  tem¬ 
perature  distributions  obtained  from  the  experiment  and  model 
in  Figs.  6  and  7,  respectively,  are  in  good  agreement  with  each 
other.  The  temperature  near  the  current-collecting  tab  of  the  pos¬ 
itive  electrode  is  higher  than  that  of  the  negative  electrode.  This 
is  attributed  to  the  electrical  conductivity  of  the  active  material  of 
the  positive  electrode  being  much  lower  than  that  of  the  negative 
electrode,  even  though  the  current  flows  near  the  tabs  of  both  the 
positive  and  the  negative  electrodes  are  similarly  high.  The  maxi- 
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Fig.  7.  Temperature  distributions  based  on  experimental  IR  image  and  modelling 
for  LIB  at  charge  times  of  (a)  3  min,  (b)  6  min  and  (c)  24  min  during  CC-CV  charging 
with  CC  charge  at  5  C. 
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Fig.  8.  Comparison  of  (a)  maximum  temperatures  and  (b)  minimum  temperatures 
obtained  from  experiment  and  model  for  LIB  during  CC-CV  charging  with  CC  charge 
at  3  and  5  C.  Only  modelling  results  shown  for  CC  charge  at  1  C. 


Fig.  9.  Comparison  of  temperature  distributions  obtained  from  experiment  along 
vertical  line  passing  through  middle  of  current-collecting  tab  from  positive  elec¬ 
trode  represented  by  arrow  in  Fig.  2  with  those  obtained  from  model  during  CC-CV 
charging  with  CC  at  (a)  3  C  and  (b)  5  C. 


mum  and  minimum  temperatures  obtained  from  the  experimental 
measurements  and  those  predicted  by  the  model,  are  shown  in 
Fig.  8(a)  and  (b)  respectively.  The  maximum  and  minimum  tem¬ 
peratures  increase  during  CC  charging  and  then  decrease  during  CV 
charging,  because  the  heat  generation  rates  decrease  rapidly  due  to 
the  rapid  decrease  of  the  charge  current  during  CV  charging.  The 
maximum  and  minimum  temperatures  obtained  from  the  experi¬ 
ment  and  modelling  are  in  good  agreement  with  each  other  over  the 
whole  range  of  DoDs  at  charge  rates  of  3  C  and  5  C  during  CC  charg¬ 
ing.  It  is  not  easy,  however,  to  resolve  experimentally  the  increasing 
and  decreasing  patterns  of  the  maximum  and  minimum  tempera¬ 
tures  during  CC-CV  charging,  when  the  rate  of  CC  charging  is  1  C, 
because  the  temperature  rise  of  the  battery  cell  compared  with 
the  atmospheric  temperature  is  not  as  clear  as  those  of  CC  charg¬ 
ing  observed  at  rates  of  3  and  5  C.  Therefore,  only  the  maximum 
and  minimum  temperatures  predicted  by  the  model  are  shown  in 
Fig.  8(a)  and  (b)  in  the  case  of  CC  charging  at  1  C.  To  compare  the 
details  of  the  temperature  distributions  obtained  from  experiment 
and  model  quantitatively,  the  temperature  distributions  obtained 
by  experiment  along  a  vertical  line  passing  through  the  middle  of 
the  current-collecting  tab  from  the  positive  electrode,  which  is  rep¬ 
resented  by  a  red  arrow  in  Fig.  2,  are  compared  with  those  obtained 
from  the  model  in  Fig.  9,  because  the  temperature  gradient  along 
that  line  is  expected  to  be  higher  than  that  in  the  other  region  of 
the  battery.  In  the  case  of  CC  charging  at  3  and  5  C,  the  temperature 
distributions  obtained  by  experiment  are  in  good  agreement  with 


those  obtained  from  the  model,  although  the  later  are  relatively 
higher  than  those  obtained  by  experimental  when  near  the  end  of 
CC  charging.  These  results  suggest  that  the  modelling  presented  in 
this  study  is  useful  to  predict  accurately  the  thermal  behaviour  dur¬ 
ing  charge  for  the  large-size  LIBs  used  in  EV  and  FIEV  applications. 
The  modelling  can  be  used  to  calculate  the  uneven  distribution  of 
heat  generation  rate  in  battery  cell  in  order  to  optimize  the  cool¬ 
ing  strategy  of  battery  packs  where  rapid  charging  is  needed.  The 
charge  model  presented  in  this  article  can  be  combined  with  the 
discharge  model  previously  reported  [19-21]  to  predict  transient 
behaviour  such  as  the  discharge/charge  voltage,  current,  state  of 
charge  and  uneven  temperature  distribution  of  LIB  in  real-life  driv¬ 
ing  cycles  for  EV  and  HEV  applications  instead  of  constant-current 
discharge/charge  conditions. 

4.  Conclusions 

A  mathematical  procedure  is  developed  to  study  the  thermal 
behaviour  of  an  LIB  during  charge.  The  two-dimensional  poten¬ 
tial  and  current  density  distribution  on  the  electrodes  of  an  LIB  are 
predicted  as  a  function  of  the  charge  time  during  CC-CV  charging 
using  the  finite  element  method.  By  comparing  the  experimental 
charge  curves  with  the  modelling  results  during  CC-CV  charging 
with  CC  charge  rates  of  1,  3,  and  5C,  the  modelling  of  the  poten¬ 
tial  and  current  density  distributions  on  the  electrodes  is  validated. 
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Then,  based  on  the  results  of  modelling  of  the  potential  and  cur¬ 
rent  density  distributions,  the  heat  generation  rate  as  a  function 
of  the  charge  time  and  the  position  on  the  electrodes  is  calculated 
to  predict  the  thermal  behaviour  of  the  LIB.  The  two-dimensional 
temperature  distributions  obtained  from  the  experiment  and  mod¬ 
elling  are  in  good  agreement  with  each  other  during  the  whole 
process  of  CC-CV  charging.  The  modelling  methodology  presented 
in  this  study  may  contribute  to  optimization  of  the  cooling  strat¬ 
egy  of  battery  packs  to  improve  the  performance  and  life  of  LIBs,  as 
well  as  to  ensuring  their  thermal  safety  in  applications  where  rapid 
charging  is  required. 
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